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Two-photon ElMl transition rates are evaluated for heliumlike ions with nuclear charges in the 
range Z = 50-94. The two-photon rates modify previously published lifetimes/transition rates of 
2^Po states. For isotopes with nuclear spin / 7^ 0, where hyperfine quenching dominates the 2^Po 
decay, two-photon contributions are significant; for example, in heliumlike ^^'^Os the two-photon 
correction is 3% of the total rate. For isotopes with / = 0, where the 2 ^Po decay is unquenched, 
the ElMl corrections are even more important reaching 60% for Z=94. Therefore, to aid in the 
interpretation of experiments on hyperfine quenching in heliumlike ions and to provide a more 
complete database for unquenched transitions, a knowledge of ElMl rates is important. 

PACS numbers: 31.10.-l-z, 31.30.Jv, 32.70.Cs, 32.80.-t 



I. INTRODUCTION 



II. THEORY 



In recent years, accurate calculations of decay rates of 
2 ■^Pq states of heliumlike ions have been performed for 
nuclear charges Z < 94, with ||, || and without [|j con- 
sideration of hyperfine quenching. Single photon decay 
of the 2 ■^Pq state to the ground state is strictly forbidden 
by angular momentum selection rules and only a weak El 
decay to the 2 state or an even weaker Ml decay to the 
2 '^Pi state is possible for isotopes with nuclear spin / = 0. 
For isotopes with 1^0, weak transitions induced by the 
hyperfine interaction also are possible. Until now, cor- 
rections to 2 ^Pq decay rates associated with two-photon 
ElMl decay have been neglected. For low-Z ions, those 
corrections are extremely tiny, scaling as Z^^; however, 
for ions such as U+^°, the ElMl decay rate is comparable 
to the single-photon rate as shown by Drake [Q. Two- 
photon ElMl rates are also expected to be significant 
for heavy two-electron ions other than U"'"^° and calcula- 
tions are needed. To date, the only ion with Z > 50 for 
which a measurement of hyperfine quenching has been 
made is heliumlike gadolinium for the isotope ^^^Gd, 
the ElMl correction to the lifetime is only 0.08%. 

We evaluate two-photon ElMl transition rates for 
2 3Po states of heliumhkc ions with Z = 50-94. The ElMl 
rate is a smooth function of Z that follows approximately 
the Z^^ law predicted by Drake [Q. Revised lifetimes 
of 2'^Po states for isotopes with / ^ 0, where hyperfine 
quenching occurs, are also given. Additionally, we give 
^Pq l^So rates for ions with Z = 50 — 94 having 7 = 0. 
Although ElMl transitions are very interesting for rela- 
tivistic theories, because of the very large negative-energy 
contributions to the transition amplitudes, direct mea- 
surements of ElMl transitions have not yet been made. 
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A. Two-photon ElMl transition rate 

The differential transition probability for the two- 
photon ElMl decay, after summation over the photon 
polarization states and integration over photon angles, 
takes the form 

dwFj = ^a^'iuliuldcj, (|M(1,2)|V |M(2,1)|') , (1) 

where the photon frequencies are related by energy con- 
servation, LOi + LU2 = Ef — Ei and the two-photon matrix 
element Af(l,2) is given by, 



2) = V (^1^1 (2)1/) 

^ ' ^ ^ E„ +uji- E, 

^ (^|Afl(2)|n) {n\E^{l)\ f) 

En + UJ2 — Ei 



(2) 



In this equation, i, /, and n designate initial, final, and 
intermediate states. Permutation of arguments in the 
matrix element M(l, 2) means that loi and UJ2 should be 
interchanged in the denominators and in the arguments 
of dipole matrix elements. Compared to the expression 
for the 2E1 transition rate given in Ref. ||^, there is the 
extra factor q;^/4, owing to differences between magnetic 
and electric dipole matrix elements. 

The matrix element M(l,2) is gauge invariant. The 
difference between the length and velocity forms of 
the electric-dipole moment in a local potential is 
{i \AEi{uj)\j) = a{E,-Ej-uj) {i j), where x(t^) is 

the gauge operator defined in Ref. and a is a constant 
of proportionality. Using this relation, the completeness 
of the basis, and the commutativity of Mi (0^2) and xi'-^i), 
we find that the length-velocity difference is for the first 
term 



{i\AEi{ui)\n) {n\M,iuj2)\f) 



Ei — En 



a{i\x{coi)Ah{cu2)\f) 
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and for the second term 
{i\Kh{uj2)\n) (n|A^i(^i)|/) 

Ei — En — LO2 



= -a{i\Nh{uo2)x{^i)\f). 



Adding the two terms, we obtain AM(1,2) = 0. The 
numerical calculations discussed below give precise seven- 
digit gauge invariance, in harmony with the theoretical 
prediction. 

After summation over magnetic substates, the two- 
photon matrix element can be further simplified, 



M(l,2) 



1 



E 



{i\\E^(uj^)\\n) (n||Afi(L^2)||/) 

Ei — En — UJl 



{i\\Nh{L02)\\n) {n\\Ei{LOi)\\f) 



(3) 



Ei — En — LO2 

Here [Ji] = 2Ji + 1 is the degeneracy of the initial state 



The quantity M(l, 2) is expressed in terms of reduced 
dipole matrix elements, which can be easily calculated. 
The units of the reduced magnetic-dipole matrix element, 
which sometimes differ in literature, are specified here by 
its relation to the transition rate 



2.69735 X 1013 Smi 

A^f = r^^ s 



A3 



9i 



where Smi = \{i ||Mi(a;)|| /)| , gi is degeneracy of the 
initial state, and A is wavelength in A. 

We consider two-photon decay of the ^Pq level in the 
single-configuration approximation, similar to Ref. [Q. 
The initial state is \i) = \ls2pi/2 [0] ) , the intermediate 
states are |n) = \lsns [1] ) (n 7^ 1) or \\snd [1] ), and the 
final state is the ground state |/) = |ls^ [0] ). Substitut- 
ing these states into (Eq. |^) we arrive at the expression 
for the two-photon amplitude used in our calculations: 



M(l,2) 



E 
E 
E 

n 

E 



(lgi/22pi/2[0] \\Ei{uJi)\\ls^,2ns^,2[l] ) {lsi/2nsy2[l] ||Ml(^2)||lg^/2[0] ) 

El - -Elsnsi/2[1] ~ ^1 
(l^l/22pi/2[0] ||^lK)||l.Si/2»rf3/2[l] ) (lgi/2n^3/2[l] || Mi (c^a) || [0] ) 

Ei — i?lsnd3/2[l] ~ ^1 

(lgl/22pi/2[0]||Ml(c^2)||l5l/2rtPl/2[l])(l5l/2npi/2[l]||^lK)||ls^/2[0]) 

El - Eianpi,2[l\ ~ ^1 

(lgl/22pi/2[0]||Afl(a;2)||lSl/2np3/2[l])(l5l/2nP3/2[l]||-gl(t^l)||lg^/2[0]) 
Ei - £^lsnp3/2[l] ~ ^2 

I 



(4) 



With the aid of the expressions for the reduced two- 
electron matrix elements given explicitly in p. 264], 
we can rewrite (Eq. ^) in terms of single-electron reduced 
matrix elements, leading to a result in agreement with 
that given by Drake ji). 



B. Line shape 

One problem of two-photon expressions is distinguish- 
ing between "pure" and cascade two-photon processes. 
The separation is needed for calculations of a total rate 
since the cascade process is already included in a single- 
photon decay rate but the pure two-photon process is 
not. (Actually, this provides a good definition for a cas- 
cade process.) Although for the off-resonance contribu- 
tion the cascade process is negligible, near the resonance 
it dominates. In |^, the subtraction of quadratic reso- 
nant terms was used to calculate the pure two-photon 
contribution. Since some details are not clear, we de- 
scribe in the appendix a subtraction procedure in more 



detail. We define a cascade process as the quadratic res- 
onant term with a fixed frequency in the numerator (this 
was not stated in Q and is important). This definition is 
reasonable because it gives a Lorentzian profile and after 
integration over frequencies of the resonant photon the 
expected single-photon rate, provided one resonant chan- 
nel dominates as is the case. Subtracting this cascade 
process from the total two-photon contribution, we can 
define the pure two-photon spectrum. Near resonance 
portion of the spectrum becomes closely antisymmetric 
and upon symmetric integration will contribute insignif- 
icantly. Resonant contributions can then be discarded 
placing grid points far away from resonant regions and 
the pure two-photon contribution can be trivially eval- 
uated. In the Appendix, we analyze the line shape of 
the pure two-photon spectrum near the strong resonance 
(there is also a second, weaker resonance but its contri- 
bution is less important). 
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III. NUMERICAL METHOD 

Since we wish to consider highly-charged ions, we base 
our calculations on the Dirac equation. To treat small 
correlation effects, of order 1/.Z, we use three starting 
potentials: Coulomb potential with Z = Ziam Coulomb 
potential with Z = Zjon — 1, and the model potential 
V{r) = —Z/r + vo{ls,r), where V{){ls,r) is the electro- 
static potential of a single Is electron. Among the three, 
the model potential takes into account correlations most 
completely, in that experimental energy levels are closely 
reproduced. Therefore, in our final compilation of two- 
photon rates, the model potential results are used. 

To carry out the sum over intermediate states, for each 
of the three potentials we use a B-spline basis set con- 
sisting of 40 positive-energy and 40 negative-energy basis 
orbitals for each angular momentum state constrained 
to a cavity of radius R = 92/Zion ao- The complete- 
ness of this basis is reflected by precise gauge invariance: 
for example, in the Coulomb potential Z = 92 for he- 
liumlike uranium, the reduced two-photon matrix ele- 
ments in the length and the velocity forms are M{1, 2)i = 
2.9395608 x lO^^ and M(l, 2)^, = 2.9395610 x lO^^. Neg- 
ative energy contributions, which almost for 2E1 transi- 
tions are small, are important in the case of ElMl tran- 
sitions. 

Since four of the terms in Eq. ^ have resonant en- 
ergy denominators and the two-photon transition rate is 
extremely sensitive to photon frequencies (~ cj^fj^f); 
replace the approximate lowest-order energies by accu- 
rate two-electron energies in these four terms: Ei 
Eils^^So), i;i.2s[i] ^ £;(ls2s35i), Ei,2p,,,[i] ^ 
E{ls2p^Po), and E^.^p.^^m ^ E{ls2p^Pi), where the 
accurate energies E{lsnl ^'^Lj) are taken from Plante 
Similar replacements were used by Drake Q in his 
calculations of the two-photon rate for 11^"+. 



IV. RESULTS 

Using the procedure described in previous section, we 
cornpute the two-photon transition rates according to the 
Eq(|l|) and after subtraction resonances integrate the dif- 
ferential rate to obtain the total two-photon rate. Since 
the behavior of the rates with Z is very smooth, we per- 
formed calculations only for eleven different ions. Table || 
summarizes our results for total transition rates. The 
three potentials, discussed above, give results in close 
agreement, but the model potential for which length- 
velocity agreement ranges 0.4-1.5% and so does approxi- 
mately accuracy should be considered the most accurate. 
(The precise gauge-invariance is achieved when denomi- 
nators are not modified in all three cases.) For U^°+our 
Coulomb Z potential result (5.567 x 10^ s~^) agrees well 
with Drake's @ (5.60 x 10^ s^^) obtained in the same 
potential. There is also small E2M2 contribution equal 
to (0.025 X 10^ s"i) which was calculated by Drake H] 



TABLE I: ElMl total transition rates for heliumlike ions in 
the range Z = 50 — 92. The rates were calculated using three 
different potentials: the Coulomb Z, the Coulomb Z—1, and a 
model potential V — —Z/r + vo{ls,r). In the model potential 
case, the difference between results obtained in length and 
velocity forms are also shown as a measure of accuracy of 
calculations. 



Z 


Coul Z 


Coul Z-1 


MP 


(L-V)/L 


50 


2.939[6] 


2.798[6] 


2.888 [6] 


1.7% 


54 


7.476 [6] 


7. 130 [6] 


7.354[6] 


1.3% 


58 


1.784 [7] 


1.706 [7] 


1.758 [7] 


1.0% 


62 


4.032[7] 


3.862[7] 


3.976[7] 


0.8% 


66 


8.686[7] 


8.335[7] 


8.575[7] 


0.7% 


70 


1.795 [8] 


1.725[8] 


1.774 [8] 


0.5% 


74 


3.574[8] 


3.442[8] 


3. 536 [8] 


0.5% 


79 


8.090[8] 


7.808 [8] 


8.014[8] 


0.5% 


82 


1.291[9] 


1.247 [9] 


1.279[9] 


0.4% 


86 


2. 355 [9] 


2. 280 [9] 


2.337[9] 


0.4% 


92 


5.567[9] 


5.403[9] 


5. 530 [9] 


0.5% 



1.6 
„ 1.5 
b 1.4 

To 1.2 

DC 

1.1 




■ Coul Z (calc) 

Coul Z (fit) 
' Coul Z-1 (calc) 

Coul Z-1 (fit) 
I Mod pot (calc) 

Mod pot (fit) 



50 60 70 80 90 
Nuclear Charge Z 

FIG. 1: Polynomial interpolations of ^-dependences for total 
transition rates calculated in the three starting potentials. 



for U^°"'". We will neglect it in all our calculations. 

In the table, we provide one extra digit beyond accu- 
racy for better interpolation. For convenience, we also 
give the interpolation polynomials that fit our calculated 
rates with sufficient accuracy. In the case of the Coulomb 
Z potential the polynomial is Z^^{1.28 - 6.4 x IQ-^Z + 
9.60x 10~^Z^); in the case of the Coulomb Z—1 potential 
the polynomial is Z^^{1.22-6AxlO-^Z + 9.8 x lO^^Z^); 
and in the case of the model potential the polynomial is 
^12(1.29 - 7.73 X lO-^Z + 1.10 x IQ-'^Z^). The rates 
in all three cases as well as the fitting polynomials are 
shown in Fig. |l|. 

The upper curve corresponds to the Coulomb Z poten- 
tial, the lower curve to the Coulomb Z—1 potential, and 
the middle one to the model potential. In addition to 
total rates we also show in (Fig. |^) the frequency distri- 
bution of the differential two-photon transition rate for 
Z = 92, Z = 79, Z = 62, and Z = 50. In the case of ura- 
nium, Z = 92, the shape resembles very closely the distri- 
bution given by Drake Q . The continuous evolution with 
Z can be seen: the maximum at the center monotonously 
decreases with Z, but the wings are rapidly developing 
for lower Z ions. 

In the case of Z = 50, the wings are large; therefore we 
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TABLE II: ElMl two-photon cor 
and ElMl rates Aeimi (1/ns) for 
decays. 



rections to lifetimes r (ps) 
ions with hyperfine induced 



TABLE III: ElMl two-photon corrections to transition rates 
A (1/s) for ions without hyperfine induced decays. 



Ion 


z 


7"old 




Aeimi 


% change 


'^''Gd 


64 


13.57 


13.56 


0.059 





08% 


^^^Tm 


69 


8.487 


8.476 


0.149 





13% 




70 


1.917 


1.916 


0.178 





03% 


173 Yb 


70 


2.113 


2.112 


0.178 





04% 


177 jj£ 


72 


1.603 


1.602 


0.253 





04% 


179 jj£ 


72 


2.556 


2.554 


0.253 





06% 


183-^ 


74 


25.54 


25.31 


0.355 





91% 


i^^Os 


76 


59.83 


58.11 


0.496 


2 


97% 


i^^Os 


76 


1.546 


1.545 


0.496 





08% 


191jj. 


77 


25.07 


24.71 


0.584 


1 


46% 


193 J J. 


77 


21.70 


21.43 


0.584 


1 


27% 


195p^ 


78 


0.9433 


0.9430 


0.686 





06% 


^^■'Au 


79 


23.04 


22.62 


0.805 


1 


86% 


199Hg 


80 


1.248 


1.247 


0.943 





12% 


201 TT . 

Hg 


80 


1.782 


1.779 


0.943 





17% 


203rpj 


81 


0.1215 


0.1210 


1.10 





01% 


205'pj 


81 


0.1192 


0.1190 


1.10 





01% 


207pb 


82 


0.8357 


0.8350 


1.29 





11% 


209Bi 


83 


0.0041 


0.0040 


1.50 





00% 




88 


5.079 


4.999 


3.14 


1 


60% 


229 rpj^ 


90 


1.955 


1.939 


4.18 





82% 


235y 


92 


2.623 


2.585 


5.54 


1 


45% 


239 


94 


3.707 


3.610 


7.29 


2 


70% 



studied the line shape more carefully with a refined scale. 
The discussion of the behavior near the origin, where res- 
onances are located, is given in the Appendix. From this 
discussion follows that the integrated contribution from 
that region is not much different from that of any other 
part of the spectrum. Since we discovered that ElMl 
corrections are significant for ions in the range Z=50-94, 
we also correct previously calculated values of lifetimes 
and transition rates for those ions. Tables |l| and [II show 
previously calculated and corrected rates. Tablc^ gives 
revised lifetimes for isotopes that have hyperfine induced 
transitions. We show only the cases where modification 
changes quoted digits. Table [II shows transition rates 



without hyperfine mixing, when corrections owing to El- 
Ml transitions are particularly large. 

The difference reaches 57.5% for Z = 94. (Drake's 



z 


Ao\d 


Aeimi 


Anew 


% change 


50 


1.363[9] 


2.87[6] 


1.366 [9] 


0.21% 


51 


1.439[9] 


3. 66 [6] 


1.443 [9] 


0.25% 


52 


1.519[9] 


4.63[6] 


1.524 [9] 


0.30% 


53 


1.605 [9] 


5. 84 [6] 


1.611[9] 


0.36% 


54 


1.695 [9] 


7.33[6] 


1.702 [9] 


0.43% 


55 


1.790 [9] 


9.16[6] 


1.799 [9] 


0.51% 


56 


1.891[9] 


1.14 [7] 


1.902 [9] 


0.60% 


57 


1.997[9] 


1.42 [7] 


2.011[9] 


0.71% 


59 


2.228[9] 


2.16[7] 


2. 250 [9] 


0.97% 


60 


2.353[9] 


2. 66 [7] 


2.380[9] 


1.13% 


62 


2.623[9] 


3.97[7] 


2.663[9] 


1.51% 


63 


2.768[9] 


4.84[7] 


2.816[9] 


1.75% 


64 


2.925[9] 


5. 87 [7] 


2. 984 [9] 


2.01% 


65 


3.090[9] 


7. 11 [7] 


3.161[9] 


2.30% 


66 


3.263[9] 


8. 59 [7] 


3. 349 [9] 


2.63% 


67 


3.446[9] 


1.03 [8] 


3. 549 [9] 


3.00% 


68 


3.638[9] 


1.24 [8] 


3. 762 [9] 


3.41% 


69 


3.839[9] 


1.49 [8] 


3.988[9] 


3.88% 


70 


4.049[9] 


1.78 [8] 


4.227[9] 


4.39% 


71 


4.273[9] 


2.12[8] 


4.485[9] 


4.97% 


72 


4.507[9] 


2. 53 [8] 


4. 760 [9] 


5.60% 


73 


4.754[9] 


3. 00 [8] 


5. 054 [9] 


6.31% 


74 


5.011[9] 


3. 55 [8] 


5. 366 [9] 


7.09% 


75 


5.287[9] 


4.20[8] 


5.707 9 

L J 


7.95% 


76 


5.574[9] 


4.96[8] 


6. 070 [9] 


8.89% 


77 


5.875[9] 


5. 84 [8] 


6.459[9] 


9.94% 


78 


6.187[9] 


6. 86 [8] 


6.873[9] 


11.09% 


79 


6.315[9] 


8.05[8] 


7.120[9] 


12.75% 


80 


6.852[9] 


9.43[8] 


7.795[9] 


13.76% 


81 


7.215[9] 


1.10[9] 


8.318[9] 


15.28% 


82 


7.663[9] 


1.29 [9] 


8. 950 [9] 


16.79% 


83 


7.971 9 


1.50 9 


9.470[9] 


18.81% 


88 


1.004[10] 3.14[9] 


1.318[10] 


31.31% 


90 1.093[10] 4.18[9] 


1.511[10] 


38.28% 


92 1.181[10] 5.54[9] 


1.735 [10] 


46.87% 


94 1.267[10] 7.29[9] 


1.996[10] 


57.51% 



prediction in uranium [Q was 46%.) Depending on ex- 
perimental situation, any ion shown in the table can be 
potentially used for the measuring ElMl transition rate. 



1.5 

S 
c 

I 0.5 

T3 



0.0 




0.2 0.4 0.6 0.1 



FIG. 2: Differential rate for two-photon decay 



CONCLUSIONS 



In this paper, we have calculated ElMl transition rates 
for heliumlike ions {Z = 50 — 94). The previous lifetime 
values are improved with adding ElMl contributions. 
The large ElMl corrections in heavy ions can be tested 
in future experiments and the ElMl transition rates can 
be extracted from experiments to study those essentially 
relativistic transitions. 
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APPENDIX: ANALYSIS OF THE LINE SHAPE 

The two-photon transition matrix element M(l, 2) has 
two poles: one occurs for n = 2 in the first term of 
Eq. (|), when wi equals the 2^Pq 2^Si (£'i) tran- 
sition energy, and the second occurs for n = 2 in the 
third term term, when wi equals the 2 ^Pq ^ 2 '^Pi (Mi ) 
transition energy. To obtain the correct line shape, the 
denominators in each of these pole terms should be mod- 
ified to include widths of the levels: En — > En — iT/2. 
(If the small width of the initial state is neglected, then 
r = r„, the width of the nth level only). The insertion 
of the imaginary term leads to an almost Lorentzian pro- 
file near the resonance. By subtracting the Lorentzians 
corresponding to the cascade processes or simply single- 
photon decays, we obtain the purely two-photon decay 
rate and its spectral distribution. After subtraction, the 
remaining contribution to the two-photon rate near the 
resonance region becomes very small when the resonant 
intermediate level is very narrow, as is the case for the in- 
termediate 2 •^S'l level and the resonances can be ignored 
in the total two-photon rate. To check that resonant 
contributions are small, we made fine grid in the reso- 
nant region and calculated the resonance contributions. 
The cascade contribution is large and is calculated sep- 
arately, using more precise methods. A precise rate for 
the 2^Pq ^2^Si {El) transition is given in Ref. §. For 
comparison with experiment, our two-photon line shape 
should be added to the Lorentzian obtained using the ac- 
curate line width and the single-photon transition rate. 
In this paper, we will show only pure two-photon profile, 
i.e., what remains after resonances are subtracted. 

The analysis of the relative contribution to the two- 
photon rate from the resonant region, using limit that 
the width is very small, gives the following expression: 



■ln|w - 1| 



1 



In this formula the unit of frequency is chosen the res- 
onant frequency ujres, li^) is the contribution obtained 
by integration from zero to to, Ii{lui) is the contribution 
with the symmetric integration from —uji to +uji with 
wi < 1. The integrals /i(wi) and I{oj) are defined be- 
low, and the expansion over small parameter F have been 
carried out: 



/i(cji) = A 




0.01 0.02 0.03 
fractional energy 



0.04 



FIG. 3: Near-resonant portion of pure two-photon spectrum 



I{uj) = A 







-4 



- u> -H —u) 
4 



(A.l) 



The illustration of the shape of the near resonance dif- 
ferential rate is given in Fig. Here the case Z — 
is considered. The differential rate near the resonance is 
very large, of order lO**, and is not shown here. However, 
owing to a precise cancellation, the total contribution af- 
ter integration is quite small. 

The symmetric integral in the near resonance region 
owing to very precise cancellation provides better nu- 
merical accuracy than asymmetric, but far away from the 
resonance when wi > 1, instead of the symmetric integral 
the asymmetric integral has to be used. We assume that 
only one strong resonance near zero where two-photon 
differential rate is proportional to uj'^ is present. The 
integral /i(a;i) is proportional to the frequency wi and 
to the width of the resonance, and the smaller the width 
the smaller resonant contribution I{u>) will be. Using the 
above equations, it can be shown that the residual reso- 
nant contribution to the two-photon rate is insignificant 
compared to the smooth continuum contribution. 

The shape in the near resonance region can be de- 
scribed approximately by the following equation 

j{x) ~ — — \- background 

[x - 1) 

if {x — 1) ^ r. The background is the smooth contin- 
uum contribution owing to all terms except the resonant 
term. There is maximum at a; = 0.012 owing to two- 
photon continuum for small x behaving like x^{\ — xY' . 
Depending on the position of most contributing levels, 
fc = 1 -f- 3 . We show the case Z = 50 for which the wings 
are the most pronounced. At the very origin, the differ- 
ential rate rapidly grows x^ if x is much less than the any 
denominator of the terms with large contributions. 

The sharp left edge of the wings is owing to the 
sharp u!^ dependence of the differential rate near the 
origin. The located very close to the origin subtracted 
resonances, which have asymmetric shapes and narrow 



6 



widths, are shown in Fig. ^ since their differential rate fast that they have no effect on the shown region as well 
near the resonance frequency is very high and needs dif- on the total rate, 
ferent scale. But resonant contributions also decrease so 
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